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THE  LIMIT  ANALYSIS  AND  DESIGN  OF  TENSION-FIELD  BE. 

By 

H.  G.  Hopkins 

Division  of  Applied  Mathematics,  Brown  University 


| 

Abstract 

| 

The  purpose  of  this  paper  is  to  develop  theoretical 
analysis  that  may  be  applied  to  discuss  the  general  prob- 
lem of  the  behavior  of  tension-field  beams.  The  analysis 
is  based  upon  the  application  of  the  theory  of  limit 
analysis  to  a plastic-rigid  model  that  approximates  the 
actual  behavior  of  such  a beam.  It  is  shown  that 
Wagner’s  tension-field  beam  analysis  is  a special  case 
of  the  present  analysis.  Apart  from  some  brief  remarks 
on  the  immediate  steps  required  in  the  development  of 
the  investigation,  no  further  discussion  of  the  problem 
is  attempted  here. 
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1.  Introduction. 

The  general  problem  of  the  determination  of  the  load- 
carrying  capacity  of  tension-field  beams  is  one  of  great  impor- 
tance  in  aircraft  structural  engineering.  These  members  are 
primary  components  in  the  wings  and  fins  of  many  types  of  air- 
craft, and  also  occur  elsewhere  in  sharply-tapered  form  as 
frames  in  the  hulls  of  flying  boats.  Much  research,  both  from 
experimental  and  theoretical  standpoints,  has  been  directed  to- 
wards the  better  understanding  of  this  problem.  The  literature 
that  describes  this  research  is  very  extensive.  In  the  following 
brief  summary,  the  references  given  are  intended,  in  the  main, 
merely  to  be  representative. 

It  should  be  noted  first  that  a tension-field  beam 
comprises  two  stiff  flanges  which  are  usually  parallel  or  nearly 
so,  a thin  web,  and  a number  of  stiffeners  which  are  usually 
evenly-spaced  and  positioned  at  right  angles  to  the  flanges. 

Such  a beavi,  regarded  as  a cantilever  with  a concentrated 
vertical  shear  load  applied  at  its  free  end,  is  illustrated  in 
Fig.  1.  The  primary  function  of  the  flanges  is  to  carry 
compressive  and  tensile  loads,  and  thereby  to  react  the  bending 
moment  at  a section.  Again  the  main  purpose  of  the  web  is  to 
carry  shear,  and  thereby  to  react  the  applied  shear  load.  Lastly, 
although  the  stiffeners  delay  the  buckling  of  the  web,  their  most 
important  function  is  to  maintain  the  separation  of  the  flanges. 
Up  to  the  time  that  the  web  buckles,  the  beam  may  be  correctly 
regarded  as  a conventional  beam  - with  shear  as  well  as  bending 
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stiffness  - and  treated  accordingly  (see,  for-  example,  [1]  )? 
However,  after  buckling  of  the  web  has  occurred  the  problem 
assumes  a radically  different  aspect.  Experiments  made  on 
typical  beam  designs  (for  which  the  flange  stiffness  is  relatively 
high)  show  that  when  the  web  buckles  it  develops  a scries  of 
waves  running  roughly  at  40°  to  the  flanges,  and  that  the  overall 
shear  stiffness  of  the  beam  drops  at  once  to  about  three-quarters 
of  its  former  value.  Prior  to  the  buckling  of  the  web,  the 
tensile  and  compressive  stresses  in  the  web  are  approximately 
equal;  but  after  buckling  of  the  web  has  occurred  the  tensile 
stress  continues  to  increase  roughly  at  the  same  rate  as  before 
whereas  the  compressive  stress  continues  to  increase  only  much 
more  slowly.  Subsequent  to  buckling  of  the  web,  the  increase  in 
applied  shear  load  is  carried  almost  entirely  by  the  tensile 
stress  in  the  web,  but  as  before  the  bending  moment  at  a section 
is  still  mainly  reacted  by  compressive  and  tensile  loads  in  the 
flanges.  Other  features  appear  with  buckling  of  the  web.  Thus 
the  flanges  become  subjected  to  distributed  transverse  loads 
from  the  web,  and  their  tendency  to  approach  one  another  under 
the  action  of  these  loads  is  resisted  by  the  stiffeners  in  which 
compressive  stresses  are  set  up.  Moreover  these  transverse  loads 
tend  to  bow  the  flanges  between  the  points  of  stiffener  attach- 
ment, and  accordingly  secondary  bending  moments,  v/hich  arc 
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greatest  near  the  points  of  stiffener  attachment,  are  induced 
in  the  flanges.  The  web  would  also  act  on  the  stiffeners  at  the 
ends  of  the  beam  in  much  the  same  way  were  it  not  for  the  fact 
that  in  practice  the  beam  design  is  almost  certain  to  embody 
special  features  near  its  ends  and  these  will  modify  considerably 
the  effects  just  described.  The  buckles  tend  to  be  smoothed  out 
by  the  flanges  and  also,  although  to  a lessor  extent,  by  the 
stiffeners.  As  the  shear  load  increases  the  amplitude  of  the 
buckles  grows  and  their  wavelength  decreases.  Ultimately  failure 
of  the  beam  must  occur,  and  this  can  happen  in  any  one,  or  a 
combination,  of  a number  of  ways  depending  upon  the  detailed 
design  of  the  beam.  For  example,  the  web  may  tear  away  from  the 
rivets  that  attach  it  to  the  flanges;  the  compression  flange  may 
fail  by  instability;  or  the  stiffeners  may  fail  by  instability 
under  the  compression  stresses  induced  in  them  by  the  web,  and 
then  follow  the  bucMes  in  the  web.  In  any  design  that  is 
satisfactory  from  a weight  point  of  view  the  various  competitive 
types  of  collapse  should  occur  together  at  failure.  This  condi- 
tion is  quite  often  approximated  through  empirical  design  modifica- 
tions following  tests.  In  any  event  it  is  important  to  note  that 
buckling  in  the  web  is  usually  very  highly-developed  indeed  prior 
to  failure  of  the  beam  - the  nominal  shear  stress  in  the  web  at 
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failure  may  be  as  much  as  twenty  times  the  buckling  shear  stress 
of  the  web,  and  approaching  the  ultimate  stress. 

The  first  attempt  to  treat  tha  post  buckling  behavior 
of  a beam  with  a thin  shear  web  is  due  to  Wagner  [ 2 ] in  his  theory 
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of  the  completely-developed  tension  field.  The  immediate  urgency 
of  this  problem  was  due  to  the  progress  in  the  design  of  aircraft, 
but  the  general  phenomena  which  occurred  had  been  noted  much 
earlier  by  Stephenson  [3]  in  experiments  on  railway  bridge 
structures.  Wagner’s  theory  assumed  that  as  soon  as  the  web  bucklei 
it  could  only  sustain  tension,  the  direction  of  this  tension 
being  of  course  parallel  to  the  main  wave  crests  and  troughs. 

Design  experience  soon  showed  that  this  hypothesis  was  too 
serious  an  over-simplication  of  the  situation  on  which  to  base 
a really  adequate  analysis.  In  particular  the  compressive  stress 
in  the  web  could  not  be  entirely  neglected.  An  appreciation  of 
this  fact  marked  the  next  stage  in  the  development  of  the  theory, 
and  gave  rise  to  the  idea  of  the  incompletely-developed  tension 
fixed.  The  new  theory  assumed  that  there  was  a compressive  stress 
in  the  web  acting  in  a direction  perpendicular  to  the  direction 
of  the  waves,  and  that  its  magnitude  was  independent  of  the 
amount  of  buckling.  Experiments  designed  to  check  the  theory 
soon  showed  however  that  the  compressive  stress  in  the  web 
continued  to  increase,  although  much  more  slowly  than  did  the 
tensile  stress,  with  the  development  of  buckling.  The  inherent 
difficulty  of  the  problem  now  led  to  the  development  of  improved 
semi -empirical  theories  (see  [4]  and  [5]).  Such  theories  have 
been  used  as  a basis  for  minimum  weight  analyses  (sec  [6  ]). 

The  task  of  developing  a complete  theory  of  tension- 
field  beams  through  an  analysis  which  treats  the  post-buckling 
behavior  of  the  web  according  to  the  von  Karmen  equations  (see  [7]) 
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for  the  large  elastic  deflections  of  thin  plates  has  proved  up 
to  the  present  time  to  be  too  difficult  a task.  This  is 
evidenced,  for  example,  by  the  attempt  of  Leggett  and  Hopkins  [8] 
to  approximate  such  an  analysis  through  an  energy  method  of 
approach  based  upon  certain  simplifying  assumptions.  The 
theoretical  results  obtained  by  these  authors  were  compared  with 
experimental  results  by  Crowther  and  Hopkins  [9]  . 

In  the  absence  of  a complete  theory  of  the  problem 
reliance  is  necessarily  mainly  to  be  made  upon  the  approximate 
and  empirical  analyses  based  upon  refinements  of  the  analysis 
initiated  by  Wagner  and  upon  empirical  conclusions  drawn  from 
experimental  observations.  This  is  not  to  say  that  the  present 
position  is  altogether  unsatisfactory  so  far  as  conventional 
designs  of  tension-field  beams  are  concerned.  However  the  lack 
of  a general  theory  does  necessarily  limit  the  understanding  of 
the  problem. 

In  the  present  paper  an  alternative  approach  is  offered 
to  the  analysis  of  the  general  problem  which  it  is  believed  will 
yield  useful  results  particularly  in  the  important  field  of 
minimum  weight  design.  Although  the  present  discussion  of  the 
application  of  the  general  theory  developed  here  is  with  refer- 
ence to  beams  with  parallel  flanges,  this  is  not  a necessary 
restriction  (see  [2  ]).  The  notable  success  of  the  methods  of 
limit  analysis  and  design,  particularly  in  their  application  to 
barred-  and  framed-structures,  has  suggested  and  encouraged  their 
use  in  the  present  study.  A brief  summary  of  the  literature 
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available  on  limit  analysis  and  design  will  now  be  given  before 
we  proceed  to  consider  the  analysis  of  the  present  problem  in 
more  detail. 

The  reader  is  referred  to  the  books  by  Van  den  Broek[10l 
and  by  Prager  and  Hodge  [11]  for  detailed  accounts  of  the  methods 
of  limit  analysis  and  design  together  with  their  various  applica- 
tions. Briefly  limit  analysis  (or  design)  in  structural  engineer- 
ing may  be  defined  as  a procedure  in  which  the  analysis  (or 
design)  of  the  structure  is  based  upon  the  ultimate  behavior  of 
the  structure.  Prager  [12]  has  defined  limit  analysis  as  being 
concerned  with  the  estimation  of  the  load  intensity  at  which  a 
given  statically-indeterminate  structure  ceases  to  be  service- 
able; and  limit  design,  on  the  other  hand,  as  being  concerned  with 
the  allocation  of  sufficient  local  yield  strength  to  the  various 
parts  of  the  structure  in  such  a manner  that  this  structure 
remains  serviceable  under  given  conditions  of  loading.  The 
basic  concepts  of  limit  design  may  be  traced  back  to  the  in- 
troduction of  ductile  materials  - wrought  iron  and  mild  steel  - 
in  engineering  structures.  However  such  concepts  were  only 
formalized  some  three  decades  ago.  The  earliest  applications 
were  to  beams  but  later  frames  were  also  successfully  treated. 
Following  the  work  on  extremum  principles  in  plasticity  by 
Drucker,  Greenberg  and  Prager  (see  [11]  ) and  by  Hill  [13] 
considerable,  and  continuing,  stimulus  has  been  given  to  applica- 
tions to  more  complicated  structural  elements  such  as  plates  and 
shells  (see  [14-]  ).  As  Hill  [13]  has  pointed  out  the  principles 
of  limit  analysis  are  most  conveniently  formulated  for 
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a perfectly-plastic  material  that  is  rigid  wherever  the  stress 
is  below  the  yield  limit. 

The  present  study  is  usefully  compared  with  the  recent 
related  work  by  Onat  and  Shield  [15]  and  by  Leth  [16]  who  discuss 
the  effect  of  shearing  forces  on  the  load-carrying  capacity  of 
wide  beams  and  of  I-beams,  respectively.  This  comparison  shows 
that  similarities  exist  in  respect  of  the  adoption  of  the  methods 
of  limit  analysis,  but,  on  the  other  hand,  dissimilarities  exist 
in  respect  of  the  correspondence  between  the  actual  physical 
structure  and  the  mathematical  model  due  to  the  present  need 
to  simulate  the  shear  buckling  of  the  web. 

The  object  of  the  present  paper  is  to  develop  a theo- 
retical analysis  that  may  be  applied  to  discuss  the  general 
problem  of  the  behavior  of  tension-field  beams.  The  theory  is 
based  upon  the  application  of  limit  analysis  to  a plastic-rigid 
model  that  approximates  the  actual  mechanical  behavior  of  such 
a beam.  The  difficulties  in  the  choice  of  a suitable  mathematical 
model  arise  chiefly  in  two  ways.  First  there  is  the  need  to 
account  for  buckling  of  the  web  and  the  relatively  small  contribu- 
tion of  the  web  compressive  stress  in  reacting  the  applied  shear 
load  once  buckling  is  well-developed.  Second  there  is  the  need 
to  account  for  strain-hardening  effects  in  the  light-alloy 
material  of  which  the  beam  is  normally  expected  to  be  made.  In 
respect  of  the  former  difficulty  it  is  of  special  interest  that 
a somewhat  analogous  situation  occurs  in  certain  problems  of 
soil  mechanics  and  in  the  design  of  masonry  structures  where  the 
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the  material  is  unable  to  take  appreciable  tension  so  that 
applied  loads  must  be  reacted  almost  entirely  through  compressive 
loads  set  up  in  the  material.  The  methods  of  limit  analysis  have 
been  successfully  applied  to  provide  solutions  of  problems  in 
these  fields  on  the  basis  of  the  assumption  that  the  material  is 
unable  to  take  any  tension  (see  [17]  and  Ll8  ] ).  The  assumption 
of  zero  tensile  yield  stress  is  not  essential  but  its  adoption 
does  greatly  simplify  the  analysis  with  comparatively  little 
loss  in  the  accuracy  of  the  values  of  the  limit  loads.  Accord- 
ingly, as  originally  proposed  by  Wagner  [2  ] , it  will  be  sup- 
posed that  the  web  is  unable  to  carry  any  compressive  stress. 

In  respect  of  the  latter  difficulty  previously  mentioned,  the 
effects  duo  to  strain-hardening  may  be  taken  into  account  through 
an  approximate,  although  reasonably  effective,  manner  by  the 
choice  of  an  effective  yield  stress  somewhat  higher  than  the 
yield  point.  The  reader  is  referred  to  a recent  paper  by 
Dwight  [19]  for  a discussion  of  a plastic  design  method  for 
aluminum  structures.  It  does  not  appear  practicable  to  assess 
directly  the  error  involved  in  making  these  simplifications  in 
the  actual  mechanical  behavior  of  the  beam,  and  the  justification 
of  results  founded  upon  the  present  analysis  must  ultimately  be 
based  upon  experimental  results. 

Finally  the  reader  should  note  that  the  relative 
simplicity  of  limit  analysis  is  due  to  the  fact  that  an  elaborate 
and  complicated  elastic-plastic  analysis  is  avoided.  Further- 
more even  when  exact  values  of  the  limit  loads  are  difficult  to 


- 


Y * 


All-114 


10 


obtain,  reasonably  close  upper  and  lover  bounds  to  these  loads 
may  be  found  through  the  use  of  the  limit  design  theorems  of 
Drucker,  Greenberg  and  Prager  (see  [11]  ).  These  theorems  will 
not  be  used  in  the  present  paper,  but  in  simple  terms  they  mean 
that  the  structure  withstands  the  applied  loads  through  the 
optimum  distribution  of  internal  stresses  and  that  the  limit 
load  is  reached  as  soon  as  a mechanism  of  collapse  is  available. 

The  general  theory  that  is  developed  here  applies  in 
slightly  modified  form  also  to  the  problem  of  the  strength  of 
composite  wall  structures  that  are  composed  of  a mild  steel 
skeleton  reinforced  by  bricks.  The  modification  is  due  of  course 
to  the  fact  that  the  brick  cannot  sustain  appreciable  tension 
although  it  will  sustain  compression.  This  application  is  not 
discussed  further  in  the  present  paper. 

2 . Notation. 

Part  of  the  notation  to  be  used  is  now  described  (see 
Fig.  2).  Let  the  web  be  rectangular  and  be  supposed  vertical 
with  one  pair  of  opposite  edges  horizontal.  Let  the  flanges  be 
rectangular  in  cross-section,  one  principal  axis  of  a cross- 
section  being  vertical  and  lying  in  the  plane  of  the  middle  sur- 
face of  the  web.  The  stiffeners  are  supposed  vertical,  at 
constant  pitch  and  symmetrically  attached  to  the  web.  Lot  0(x,y) 
be  a rectangular  Cartesian  frame  of  reference  in  the  middle  surface 
of  the  wob,  the  x-  and  y-axes  being  horizontal  and  vertical, 
respectively-  Let 

<L,d,t  = web  length,  depth  and  thickness 5 
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2a, b = depth  and  width  of  flange  cross-section; 
d0  = tensile  yield  stress; 

F,Q,M  = axial  and  shear  forces,  and  bending  moment  in  flange; 

Up,Vp  = velocity  components  in  neutral  axis  of  flange; 

d . d » t = stresses  in  web  averaged  across  thickness; 
x’  y’  xy 

vx,v  = velocity  components  in  web  averaged  across  thickness;  and 

y 

ex’ey’YXy  = s^ra:'-n  components  in  web  averaged  across  thickness. 


Other  notation  is  defined  when  first  introduced. 


3*  Preliminary  Analysis. 

It  is  convenient,  before  proceeding  to  the  detailed 
analysis  of  the  flanges  and  the  web,  to  describe  the  general 
approach  which  follows  closely  that  developed  by  Prager  (see  [20] 
and  [2l]  ). 

The  mechanical  behavior  of  a porfectly-plastic  rigid 
continuum  is  specified  in  terms  of  independent  generalized  strain- 
rates  e±  and  stresses  <3±( i = 1,2,..., I),  these  quantities  being 
associated  in  such  a way  that  the  rate  of  plastic  work  per  unit 
volume  is 
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The  generalized  strain-rates  arc  derived  from  generalized  veloci- 
ties v^  (j  = 1,2, ... ,J)  which  specify  the  typo  of  plastic  flow 
envisaged  and  reflect  the  basic  kinematical  assumptions  of  the 
theory.  In  general  once  the  v.'s  are  assigned  and  the  e*s  have 
been  determined,  the  q^s  arc  best  found  from  considerations  of 
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virtual  work.  The  c^'s  may  always  be  modified  in  such  a way  to 
have  the  same  dimensions  as,  for  example,  conventional  strain- 
rates.  The  same  applies  to  the  d^'s.  Indeed  if  is  modified 
to  e^/c^,  where  c^  is  some  constant,  then  d^  must  bo  modified  to 
cidi  » if  W is  to  remain  unchanged.  This  modification  is  essen- 
tial if,  as  is  often  convenient,  the  e^'s  and  d^'s  are  to  be 
regarded  as  the  rectangular  Cartesian  co-ordinates  of  points  in 
I-dimensional  strain-rate  and  stress  spaces.  Further  the  points 
e^  and  d^  , and  hence  the  states  of  strain-rate  and  stress,  may 
be  associated  with  the  vectors,  say  e and  c respectively,  having 
their  initial  points  at  the  origin  and  their  terminal  points  at 
and  d..  . Moreover  it  may  be  convenient  to  superimpose  the 
strain-rate  and  stress  spaces  in  which  case  all  e^'s  and  d^'s 
must  be  chosen  so  as  to  have  the  same  dimensions.  In  the 
absence  of  effects  such  as  strain-hardening,  viscosity  and  inertia, 
no  fundamental  time  is  involved  in  the  problem,  and  accordingly 
the  velocities  are  only  determinate  to  within  an  arbitrary  con- 
stant factor  of  proportionality. 

The  physical  laws  involved  in  the  problem  are  those  of 
equilibrium;  plastic  flow  without  fracture  or,  within  a rigid 
(elastic) region,  strain  compatibility;  and  mass  conservation. 

These  laws  are  just  sufficient  to  provide  the  necessary  set  of 
mathematical  relations  on  the  unknown  field  quantities  and  either 
the  external  forces  or  the  geometry  according  as  the  problem  is 
one  of  analysis  or  design.  The  physical  laws  provide  not  only 
field  equations  and  boundary  conditions  but  also  certain  funda- 
mental requirements  in  respect  of  continuity  on  field  quantities. 
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The  equations  of  equilibrium  involve  both  internal 
and  external  stresses.  Any  internal  stress  which  does  not 
explicitly  occur  in  W may  be  associated  with  a strain-rate  that 
is  identically  zero.  Such  a stress  is  called  a reaction.  The 
actual  derivation  of  the  equations  of  equilibrium  is  standard. 

The  yield  condition  is  expressed  by 


(3.2) 


where  the  yield  function  f assumes  negative  values  for  stress 
states  below  the  yield  limit.  In  the  simplest  case  f is  every- 
where continuously  differentiable  with  respect  to  the  d^'s,  i.e* 
the  yield  surface  has  everywhere  a continuously-turning  tangent 
plane,  and  f is  called  regular.  Otherwise  the  yield  surface  is 
called  singular,  and  points  where  all  the-^'s  do  not  exist  are 


called  singular  points.  Note  that,  although  it  is  convenient 
to  use  geometrical  language,  expressions  such  as  'yield  surface' 
must  be  interpreted  appropriately  in  any  particular  context. 

The  most  important  example  of  a singular  yield  surface  occurs 
when  the  (closed)  yield  surface  S comprises  a finite  number  of 
(open)  regular  surfaces  Sp(r  = 1,2,..,,R)  so  that  singularities 
may  occur  only  at  the  intersections  of  the  Sr's.  This  situation 
arises  when  more  than  one  physical  mechanism  is  available  to 
admit  plastic  flow  at  certain  stress  states.  In  such  a case  the 
yield  limit  Eq.(3«2)  is  expressed  in  terms  of  a finite  number  of 


me* 


7* 
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continuously  differentiable  yield  functions 

f p ( 0-, , ^2  i * * * i ) (r  — l,2,«e-,R)  (3*3) 

where  all  the  fr's  assume  negative  values  for  stress  states 
below  the  yield  limit,  and  at  least  one  f is  zero  with  no  fr 
positive  for  stress  states  at  the  yield  limit.  Except  perhaps 
at  the  boundaries  of  the  Sr's  there  is  always  a v/ell -defined 
outwards -drawn  unit  normal  ^(n^)  to  S where 

n oc  (i  = 1,2,. ..,1);  (3.4) 

i oo^ 

and,  at  a singular  point  of  S,  n is  defined  by 

af„ 

ni  °c  S ar  ^ (i  = 1,2,..., I),  (3.5) 

where  the  ar's  are  non-negative  and  not  all  zero,  with  the  sum- 
mation extending  over  all  values  of  r for  which  Sr  passes  through 
the  singular  point  in  question.  Relations  (3.4)  and  (3.5) 
coincide  at  regular  points,  and  the  latter  are  taken  as  an 
extended  definition  of  n that  applies  at  all  points. 

The  general  problem  of  the  determination  of  the  actual 
form  of  the  yield  function  is  not  discussed  here,  but  in  many 
important  cases  of  practical  interest  the  procedure  is  relatively 
straightforward.  A fundamental  requirement  is  that  the  yield 
surface  is  everywhere  concave  towards  the  origin.  However  once 
f is  known  the  formulation  of  the  flow  rule  on  the  basis  of  the 
theory  of  the  plastic  potential  due  to  von  Mises  is  immediate. 
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This  theory  demands  that,  for  plastic  stress  states,  the  plastic 


strain-rate  vector  and  the  outwards -drawn  normal  vector  to  the 


yield  surface  at  this  stress  state  have  the  sane  direction. 


Otherwise  the  plastic  strain-rate  vector  is  identically  zero. 


Thus 


e = 0 wherever  either  f < 0 or  f = 0 and  f'  < 0, 


e - \n  wherever  f = f'  = 0 (X  > 0), 


(3.6) 


\ being  a function  of  position  indeterminate  to  the  extent  of  a 


constant  multiplicative  factor.  It  is  unnecessary  to  give 


explicit  attention  to  the  conditions  on  f1  (the  prime  denoting 


time  differentiation)  because  these  are  automatically  satisfied 


as  the  incipient  stress  field  at  collapse  is  supposed  independent 


of  time  at  least  to  a first  approximation.  Equivalently,  the 
relations  (3.6)  are  expressed  through  the  relations  (3.*+)  and 


(3.5)  in  the  forms 


Cj  = 0 wherever  either  f < 0 or  f = 0 and  f1  < 0, 

> (i=l,2,...,I), 

ei  = X — wherever  f = f = 0 (X  > 0),  (3./) 

QcJi  J 


for  a regular  yield  surface,  and 


= 0 wherever  either  f < 0 or  f = 0 and  f ’ <0  (r=l,2, . . , ,R) , 


X 2 a — £.  wherever  f =f'  =0  for  at  least  one  r 
r r ad.  r r 


( X>  0;  a„  = 0 if  either  f < 0 or  f =0  and  f ' <0,  and 

r T T*  T*  • 


a > 0 if  f = r = 0), 
r - r r 


(1*1,2, ..,l), (3.8) 


\ r*  - 4-h - ■ 
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quite  generally  for  a singular  yield  surface.  In  relations  (3*8) 
the  ar’s  are  otherwise  completely  arbitrary,  and  variations  in 
their  ratios  correspond  to  different  combinations  of  plastic 
flow  mechanisms. 

Plastic-rigid  material  is  incompressible,  and  hence  the 
density  is  constant.  It  follows  that  the  physical  law  of  mass 
conservation  is  expressed  in  terms  of  a continuity  condition 
involving  only  velocity  components.  The  present  theory  is  not 
applicable  to  problems  of  fracture.  For  our  purposes  fracture 
may  be  defined  as  severe  local  plastic  deformation  resulting  in 
the  formation  of  new  surfaces  to  the  body.  Caro  is  needed  to 
distinguish  between  deformations  which  are  properly  regarded  as 
being  within  the  scope  of  plasticity  theory  and  deformations 
involving  fracture  which  are  strictly  outside  the  scope  of  this 
theory.  Even  so  plasticity  theory  is  commonly  applied  to  the 
analysis  of  problems  which  involve  either  highly-localized  regions 
of  fracture  or  incipient  fracture,  and  this  practice  extends 
considerably  the  application  of  the  theory. 

Finally  it  is  necessary  to  introduce  certain  definitions 
in  respect  of  discontinuities  of  field  quantities.  A surface  T 
of  isolated  discontinuity  is  said  to  be  of  order  n(>  1)  with 
respect  to  the  field  quantity  , if  ^together  with  all  its 


space  derivatives  up  to  and  including  those  of  order  n-1,  but 
not  all  derivatives  of  order  n,  are  continuous  across  T.  If  the 
discontinuity  effects  3~ itself  then  it  is  said  to  be  of  order 
zero.  The  basic  field  quantities  are  of  course  Vj,^  and  X. 
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On  each  side  of  T the  various  field  quantities  are  related  by- 
equations,  either  differential  or  finite  in  nature,  valid  at  all 
points.  For  example,  there  are  the  differential  equations  of 
equilibrium  and  the  yield  condition  all  involving  the  generalized 
stresses.  In  what  follows  it  is  to  be  understood  that  all  equa- 
tions are  fundamental,  and  not  derived,  in  the  sense  that  they 
express  the  various  physical  laws  in  the  simplest  possible  manner. 
Then  with  respect  to  T and  each  3 , a weak  discontinuity  has  an 
order  not  lower  than  the  order  of  the  highest  derivative  of 
occurring  in  the  physical  equations.  Otherwise  the  discontinuity 
is  called  strong,  and  is  a contact  discontinuity.  If  the  material 
in  the  immediate  neighborhood  remains  rigid  then  there  is  said 
to  bo  localized  plastic  flow.  The  above  definitions  need  to  be 
extended  when  the  time  enters  explicitly  into  the  problem. 

Further  any  (virtual)  discontinuity  is  called  natural  or  artificial 
according  as  the  physical  equations  do  or  do  not  remain  valid  in 
its  immediate  neighborhood.  Note  that  if  there  is  a dis- 
continuity in  3"  across  T then  it  is  properly  to  be  regarded  as 
the  limit  of  a continuous  distribution  of  t which  changes  by  a 
fixed  amount  across  a narrow  region  enclosing  T as  this  region 
everywhere  shrinks  up  to  T. 

4,  Analysis  of  the  Flanges. 

The  flanges  are  treated  according  to  a simple  engineer- 
ing theory  of  perf ectly-plastic  rigid  beams,  i.o.  the  flanges 
are  approximated  by  one-dimensional  continua.  The  relation  be- 
tween the  present  approximate  theory  and  the  exact  theory  is 
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reasonably  conjectured  to  be  much  the  same  in  plasticity  as  it 
is  in  elasticity.  Exactly  as  in  elastic  beam  theory,  transverse 
shearing  deformations  being  neglected,  the  fundamental  kine- 
matical  assumption  is  that  a cross-section  of  a flange  initially 
normal  to  the  neutral  axis  remains  normal  to  the  deformed  neutral 
axis.  If  Up  and  Vp  denote  the  axial  and  transverse  velocity 
components  of  particles  lying  on  the  neutral  axis,  then  the  rate 
of  extension  e and  of  curvature  x of  the  neutral  axis  are  given 


dUp 


(4.1) 


Let  F,  Q and  M denote  axial  force,  vertical  shear  force  and 
bending  moment.  Then  the  rate  of  plastic  work  per  unit  length 
of  flange  is 


W = F e + M x . 


(4.2) 


Hence  F,  M,  e , x,  Up  and  are  corresponding  generalized  stress, 
strain-rates  and  velocities,  and  Q is  a reaction. 

It  is  simple  to  establish  certain  general  conditions 
of  continuity  on  the  fundamental  field  quantities  F,  Q,  M,  Up 
and  Vp. The  equilibrium  of  an  elementary  slice  of  a flange 
enclosing  an  arbitrary  fixed  cross-section  shows,  when  the 
thickness  of  this  slice  is  made  to  tend  to  zero,  that  F,  Q,  and 
M are  continuous.  The  law  of  mass  conservation  applied  to  the 
same  slice,  remembering  that  the  density  is  con&ant  and  that 
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finite  changes  in  the  flange  cross-sectional  area  would  involve 

fracture,  establishes  that  u^  is  continuous.  A discontinuity  in 

v_  when  properly  regarded  would  correspond  to  large  transverse 
r 

shear  strain-rates  within  the  elementary  slice.  Inasmuch  as  the 
effect  of  transverse  shear  strain-rate  is  completely  ignored  in 
the  present  theory,  it  is  accordingly  necessary  for  consistency 
to  assume  that  v^,  is  continuous.  Thus  at  all  cross-sections, 


[F,Q,M,Up,Vp  ] - 0, 


(4.3) 


the  square  brackets  denoting  discontinuity  of  the  enclosed 
quantity.  The  sign  conventions  for  the  quantities  introduced 
above  are  shown  in  Fig.  2. 

The  forces  acting  on  an  elementary  slice  (of  length 
6x)  of  the  upper  flange  are  shown  in  Fig.  2,  It  follows  that 
the  differential  equations  of  equilibrium  are 


a§  - t(‘V)y-fcl  “ °’ 

to  + t(Vy=*i  = °’ 


dh 

dx 


Q 


= o, 


(4.10 

(4.5) 

(4.6) 


the  upper  flange  lying  along  the  line  y = id,  Similar  equations 
apply  to  the  lower  flange.  Here,  for  simplicity  and  with  little 
resulting  loss  in  accuracy,  the  offset  of  the  web  edge  from  the 
flange  neutral  axis  is  neglected;  this  involves  the  omission  of 
a term  - at(TXy.)y._^d  from  the  left-hand-side  of  Eq.(4.6).  If 
there  is  a distributed  downwards- acting  load  p per  unit  length 
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of  flange,  then  a term  +p  must  be  added  to  the  left-hand -side 
of  Eq.  (4,5). 

As  established  by  Onat  and  Prager  [22  ],  the  critical 
combination  of  F and  M admits  full  plasticity  at  a cross-section 
corresponds  to  stress  distributions  over  a cross-section  in  which 
axial  fibers  are  at  the  yield  stress  either  in  tension  (+dQ)  or 
in  compression  (-dQ).  It  is  straightforward  to  determine  the 
form  of  the  yield  functicn(s)  for  any  shape  of  cross -section. 

For  example,  in  the  case  of  a rectangular  shape  to  which  attention 
is  confined  here,  the  yield  limit  is  expressed  in  terms  of  two 
yield  functions,  viz. 


where 


fx  = F2/Fq2  + M/M0  - 1, 
f2  = F2/F  2 - M/M0  - 1, 


= 2abdc  , M = a"bdo 


(4.7) 


> 


Here  FQ  and  M0  are  the  limit  axial  force  and  moment  in  pure 
extension  and  bending,  respectively.  As  will  be  seen  in  a 
moment  the  converse  statement  is  not  necessarily  true.  In 
proceeding  it  is  necessary  to  modify  the  generalized  stresses 
to  F/F0  and  M/MQ  , and  the  generalized  strain-rates  to  Fc e 
and  M x . The  corresponding  generalized  velocities  can  be  taken 

dVr, 

as  F0Up  and  M0  -=— £«  (assuming  that  F0  and  M0  are  constant)  so 


that  F0 e = ^ (F0Up)  and  M0  x = iL  (mq  ^4L),  but  these  quantities 
are  not  explicitly  required  in  this  form.  Then  in  the  correspond- 
ing stress  plane  the  yield  locus  is  represented  by  two  parabolic 
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arcs  (ABC  and  CDA  in  Fig.  3)  with  singular  ooints  only  at  their 
common  points  (A  and  C).  The  flow  rule  is  given  by  the  follow- 
ing relations. 


■tegime 

Fo 

Mo 

Conditions 

A 

2 \ 

-(l-2q) X 

F/F0=l,  M=0;  0 < q < 1 

ABC 

2\F/Fc 

\ 

- 1 < F/Fq  < 1 , M >0 

C 

-2\ 

(l-2q)\ 

F/Fq=-1,  M=0;  0 ^ q < 1 

CDA 

2\F/Fq 

-1  < F/F0  <1,  M < 0 

(4.8) 


The  quantity  \(x)  is  non-negative,  and  variations  in  q(x) 
correspond  to  different  combinations  of  the  admissible  plastic 
flow  mechanisms  at  a singular  point. 

The  law  of  mass  conservation  does  not  yield  any  further 
information  that  is  directly  relevant  to  the  solution  of  the 
problem.  It  determines  only  the  rate  of  change  A'  of  cross- 
sectional  area, 

A'  = - Ae  . (4.9) 

The  solution  of  the  problem,  so  far  as  the  flanges 
are  concerned,  must  satisfy  the  continuity  solutions  (4.3),  the 
equilibrium  equations  (4.4)-(4*6)  and  the  flow  rule  (4.8).  Note 
that  the  matching  of  the  solution  obtaining  in  either  rigid  or 
plastic  regimes,  i.e.  at  rigid-plastic  or  plastic-plastic  inter- 
faces is  achieved  through  the  continuity  relations.  The  solu- 
tion, if  it  is  to  be  strictly  within  the  present  framework,  will 
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not  admit  discontinuities  in  F,Q,M,up  and  vp. 

Finally  it  remains  to  discuss  the  general  question 

of  discontinuities  in  the  field  quantities.  The  highest 

dv„ 

derivatives  of  F.QjMjUp  and  that  occur  in  the  physical 

equations  are  of  the  first  order.  A weak  discontinuity  in  these 

dF 

quantities  is  now  considered.  First  a discontinuity  in  -j™  and 
H involves  discontinuities  in  (Txy)y_^  and  (<Jy)  » 

respectively,  and  must  therefore  be  discussed  later  when  the 
analysis  of  the  web  has  been  completed  (see  p.46).  A dis- 
continuity in  ~ involves  a discontinuity  in  Q which  violates 
one  of  conditions  (4.3),  and  is  not  admissible.  Secondly  dis- 
continuities in  e and  x can  occur  only  either  within  or  at  the 
boundary  of  a plastic  regime.  Then,  remembering  that  F is 
continuous,  the  following  results  are  found  from  the  flow  rule. 


Regime 


2 (F/F  ) [\J 
o 

-2 [\]  “ 
2(F/Fo)[X] 


- LCL-2qTn 

ixT 

Kl-2q)\f 

-m 


(4.10) 


Here  [e]  and  [x]  are  supposed  not  both  zero,  and  the  results  for 
the  three  distinct  cases  follow  immediately  and  are  tabulated 
below. 
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[e] 

[h] 

Regime 

Conditions 

: . . 

^0 

ABC , CDA 
A,  C 

F,[\]  *0 

> « 
F,[\]  / 0 , [(l-2q)\ 

/o 

0 

A,  C 

[\]  / 0 , [(l-2q)\]=  0 

0 

/ 

^0 

B,D 

» 

A,C 

[X]  *0 

[\]  = 0,  [l-2q]  /O 

(4.11) 


Thus  If  either  regime  ABC  or  CDA  applies  then  discontinuities  in 
axial  strain-rate  and  curvature-rate  must  occur  together  and  are 
not  independent,  save  when  there  is  zero  axial  force  in  which 
case  no  discontinuity  is  permissible  in  axial  strain-rate  al- 
though a discontinuity  is  permissible  in  curvature-rate;  and  if 
either  regime  A or  C applies  then  discontinuities  in  axial  strain 
rate  and  curvature-rate  may  occur  separately  or  together.  At 
first  sight  the  behavior  in  regime  A or  C may  appear  rather 
surprising,  but  it  is  simply  explained  by  the  fact  that  if  all 
axial  fibers  at  a section  are  stressed  to  the  tensile  or  compres- 
sive yield  stress  it  does  not  follow  that  there  is  constant  axial 
velocity  across  the  section  - the  distribution  of  axial  ■velocity 
is  only  required  to  bo  a continuous  linear  function,  either 
uniformly  positive  or  negative,  of  the  distance  from  the  neutral 
axis,  and  hence  bending  as  well  as  extension  may  be  accommodated. 
This  completes  the  discussion  of  weak  discontinuities. 

Next  consider  strong  discontinuities.  A discontinuity 
in  F,Q  or  M violates  equilibrium,  and  is  therefore  artificial. 
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Such  discontinuities  are  not  considered.  Discontinuities  in 


d\’ 


Up  and  cannot  occur  within  a rigid  regime  but  may  occur 
elsewhere.  Here  [Up  ] and  [^—]  are  not  both  zero,  and  hence 
one  of  e and  x becomes  very  large  within  the  elementary  slice. 

Now  F remains  continuous,  and  hence  X becomes  very  large  within 
the  elementary  slice.  Let 

I = lim  J \(x)dx  (>0),  J = iim  J ^jl-^qCx)^ \(x)dx  (-1  £ J + I), 

(4.12) 


the  limits  being  taken  as  the  thickness  of  the  elementary  slice 
is  made  to  tend  to  zero.  The  following  relations  now  follow 
through  integration  of  the  flow  rule. 


Regime 

i — i 

o 

*5 

i— j 

A 

21 

-J 

ABC 

2 (F/Fq)  I 

I 

C 

-21 

J 

CDA 

2 (F/F  Q)  I 

-I 

(4.13) 


The  results  for  the  three  distinct  cases  now  follow  immediately. 
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Regime 

Conditions 

AO 

A o 

ABC, CDA 

F A o 

A,  C 

J £ 0 

AO 

0 

A,C 

J = 0 

n 

i° 

L 

A o 

B,D 

MM 

(4.14) 


Thus  if  either  regime  ABC  or  CDA  applies  then  discontinuities  in 
axial  velocity  and  angular  velocity  must  occur  together  and  are 
not  independent,  save  when  there  is  zero  axial  force  in  which 
case  no  discontinuity  is  permissible  in  axial  velocity  although 
a discontinuity  is  permissible  in  angular  velocity,  and  if  either 
regime  A or  C applies  then  a discontinuity  in  axial  velocity 
always  occurs  and  may  also  be  accompanied  by  an  independent  dis- 
continuity in  angular  velocity.  These  discontinuities  are 

classified  in  the  following  way.  First  if  [FQUp  ] = 0 and 
dv-n> 

[ Mq  0 then  the  discontinuity  is  immediately  recognized 

as  corresponding  to  a 'yield  hinge'  long  familiar  in  the  plastic 

analysis  of  beams.  Second  if  [F0Up  ] / 0 and  [M0  A 0 then 

the  discontinuity  is  called  an  'extensible  yield  hinge',  and 

was  first  introduced  by  Onat  and  Prager  [22]  . Last  if 

[FQUp  ] A 0 and  [MQ  ] =0  then  the  discontinuity  is  simply 

dx 

called  an  'extension'.  A discontinuity  in  axial  velocity  in- 
volves very  high  rates-of -change  of  cross-sectional  area  for 

i lim  J A'  dx  s - [up]  , (*+.15) 

the  limit  being  taken  in  the  usual  way.  This  situation  corres- 
ponds to  incipient  fracture  simultaneously  with  plastic  flow, 
and  is  accordingly  not  strictly  admissible  within  the  framework 
of  the  present  theory.  Moreover,  if  discontinuities  occur  in 
either  axial  velocity  or  angular  velocity  then  high  rates  of 
strain  of  axial  fibers  must  occur.  The  present  theory  does  not 
take  account  of  such  effects.  Accordingly  it  must  be  expected 
that  viscous  effects  will  immediately  serve  to  modify  such  sharp 
discontinuities.  Moreover  effects  due  to  inertia  and  strain- 
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hardening  will  soon  become  important  at  such  a section.  In  a 
general  way  the  discontinuities  will  be  spread  out  rapidly  over 
a small  length  of  the  flange.  The  discontinuities  are  artificial, 
but  nevertheless  analysis  within  the  present  framework  may 
reasonably  be  expected  to  furnish  a sufficiently  adequate 
description  of  the  actual  mechanical  behavior  except  in  the 
immediate  neighborhood  of  such  discontinuities. 


5.  Analysis  of  the  Web. 

The  web  is  treated  according  to  an  approximate  theory 

of  plane  stress  for  a hypothetical  material  which  is  unable  to 

support  compressive  stresses.  Thus  the  web  is  approximated  by 

a two-dimensional  continuum.  So  far  as  the  assumptions  of 

plane  stress  are  concerned,  the  relation  between  the  present 

approximate  theory  and  the  exact  theory  is  reasonably  conjectured 

to  be  much  the  same  in  plasticity  as  it  is  in  elasticity.  The 

generalized  stresses  are  the  mean  (i.e.,  thickness  averaged) 

stresses  a , d , x , and  the  corresponding  generalized  strain- 

* y xy 


rates  are  the  mean  strain-rates  e. 


Y-* 


'x’  y7  'xy 

denote  the  components  of  mean  velocity  then 


If  v and  v 

y 


^2 

0K 


> e 


^ * YXy  ax  Qy 


(5.1) 


The  forces  acting  on  a small  element  6 x x 6v  x t are 
shown  in  Fig.  2.  The  differential  equations  of  equilibrium  are 
therefore 


8d. 

ax" 


x 


= 0 


d%. 

ax 


xy 


9dy 

+ = 0 . 
ay 


(5.2) 


The  rate  of  plastic  work  per  unit  area  of  middle 


surface  is 


w = + tivev  + TxvYvJ- 


(5,3) 


x~x.  y y xy  * xy 

The  principal  stresses  are  denoted  by  d and  d2,  and 
are  distinguished  through  the  choice 


(5A) 


In  the  case  of  equality  in  condition  (5.4-)  there  is  stress 
isotropy.  If  cp  is  the  anti-clockwise  rotation  of  the  d^  - and 
d2  - directions  from  the  x - and  y - axes,  respectively,  then 


2<*i  = dy.  + + r , 2d 2 = d x + “ r > 

cos  2 <p  = (dx-dy)/r  , sin  2q>  = 2TXj/r  » 
r = + l (d„-d„)2  + 4x2wV  > 


> 


(5.5) 


\vY°y;  T^xyf 

(cp  is  of  course  only  determinate  to  within  a multiple  of  ft), 
and  conversely 


2 2 

dx  = cos  tpd-^  + sin  9 d2  , 


d = sin  9 d,  + cos®9  d2  > 

y 

= £ sin  2 9 (d1  - d2). 


(5.6) 


The  yield  condition  is  required  to  express  the  fact 
that  the  material  will  not  sustain  pure  conpresion  and  may  yield 
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when  either  principal  stress  achieves  the  yield  stress  dD  in 
simple  tension.  Hence  the  yield  condition  restricts  d^  and  d2 
to  lie  either  between  or  at  the  values  0 and  dQ  , i.e.,  v/ithin 
the  square  ABCD  (see  Fig,  4)  drawn  in  a stress  plane  in  which  d^ 
and  d2  are  taken  as  rectangular  Cartesian  co-ordinates.  However, 
remembering  the  condition  (5.4),  only  the  triangle  ABC  is 
strictly  relevant. 

The  yield  curve  ABCD  is  imagined  folded  about  the 
diagonal  AC,  and  consists  of  the  lines  AB  and  BC  each  described 
twice.  The  diagonal  AC,  its  end  points  being  excluded,  is  not 
part  of  the  yield  curve.  Accordingly  there  are  two  regular 
plastic  regimes  AE  and  BC,  and  three  singular  plastic  regimes 
A,B  and  C.  Thus  the  yield  condition  is  expressed  in  terms  of 
two  yield  functions,  viz. 


fl  = " d2»  f2  = dl  " 6o  ’ 


(5.7) 


subject  of  course  to  condition  (5.4).  This  yield  condition  is 
also  simply  expressed  in  terms  of  the  stresses  dx,  d and 
as  follows.  Now  for  the  regimes  A,  AB  and  B,  0 <;  d^  < dQ  and 
c?2  = 0 ; and  for  the  regimes  B,  BC  and  C,  d^  = dQ  and 
0 1 d2  - d0  * ^en  i'rom  Eqs.(5'.5)  and  (5.7), 

for  regimes  A,  A3,  B,  ^ 

(dQ-dx)  + (do”dy^  for  reGimes  B>  BC>  c» 


■f 


V = d + d 
1 x y 


(5.8) 


and  hence 


where 


3 = 


0 for  regimes  A,  AB, 

1 for  regimes  B,  BC, 


w 


(5.9) 


Note  that  regime  E is  characterised  by  either  value  of  j. 

If  f denote  the  yield  condition  then  the  flow  rule  is 
expressed  by 


+ 


)/-&— 

9Txy 


- X(Xjy)  2.  0, 


(5.10) 


the  first  order  partial  differential  coefficients  of  f being 
given  the  usual  freedom  of  interpretation  at  singular  points. 

The  law  of  mass  conservation  determines  the  rate-of- 
charige  of  the  we b thickness,  or,  equivalently,  the  mean  trans- 
verse strain-rate  e , in  terms  of  e and  el 

t*  = tez  = - t(ex  + e ) . (5.11) 

The  set  of  five  equations  (5.2),  (5.9)  and  (5.10)  in 
the  five  unknowns  <s„  . 4 , x , v„  and  v form  the  basis  for 
the  determination  of  the  stress  and  velocity  fields  in  the  plastic 
region.  Note  that  the  first  three  of  these  equations  involve 
only  the  stresses.  The  other  equations  involve  the  stresses  and 
the  velocities,  and  are  homogeneous  in  the  latter. 

In  the  rigid  (elastic)  region  the  equations  of 
equilibrium  remain  the  same,  but  now  the  strain-rates  all  vanish 
identically  and  the  yield  condition  is  replaced  by  the  requirement 
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oi  compaxxDxxxLy, 


^2  o 

(— 2 + — ^)(dv  + d,„)  = 0. 

6x  6y 


'x 


(5.12) 


The  set  of  three  equations  (5.2)  ana  (5.12)  in  the  three  unknowns 
dx  ’ an(*  Txy  f°rm  the  kas*s  for  the  determination  of  the 
stress  field  in  the  rigid  ( elastic)  region. 

It  is  necessary  to  establish  certain  conditions  on  the 
field  quantities  at  the  plastic-elastic  interface.  Let  the 
interface  be  a simple  open  curve  V with  everywhere  a 
continuously-turning  tangent  except  perhaps  at  a finite  number  of 
isolated  points.  Let  the  normal  and  tangential  directions  to  T 
be  denoted  by  increasing  values  of  local  rectangular  Cartesian 
co-ordinates  n and  3,  respectively.  The  equilibrium  of  a small 
rectangular-shaped  volume  6s  x 6n  x t shows,  as  the  thickness  6n 
is  made  to  tend  to  zero,  that  the  normal  stress  dn  and  the  shear 
stress  Tnsare  continuous,  but  that  the  tangential  stress  need 
not  be  continuous  (see  Fig.  5).  The  law  of  mass  conservation 
applied  to  the  same  volume,  remembering  that  the  density  is 
constant  and  that  finite  changes  in  the  web  thickness  would 
involve  fracture  which  is  here  excluded,  establishes  that  vn  is 
continuous  but  that  vg  may  be  discontinuous.  Thus 


tOn-'W  Tn]“  0 


(5.13) 


across  a plastic-elastic  interface,  but  d and  v may  be  dis- 

s s 

continuous  across  such  an  interface. 


■V  f 


It  will  be  established  later  that  the  directions  of 
principal  stress  at  any  point  of  a plastic  stress  field  are 
mathematically-preferred  directions.  This  fact  at  once  focuses 
attention  on  the  field  of  principal  stress  trajectories  as  the 
fundamental  unknown  element  to  be  determined  in  the  solution  of 
the  problem.  Let  the  two  orthogonal  families  of  lines  of 
principal  stress  be  identified  through  the  orthogonal  curvilinear 
co-ordinates  a and  p.  Hence  a(x,y)  and  p(x,y)  are  supposed 


continuously  differentiable  functions  such  that 

9a  d§  da  38  n 

» si  * dp  = °' 


(5.1V) 


A member  of  one  family  of  principal  stress  trajectories  is  then 
defined  by  an  equation  {3  = const,  and  is  called  an  a-line;  and 
similarly  a member  of  the  other  family  is  defined  by  an  equation 
a = const,  and  is  called  a p-line.  The  two  families  are  dis- 
tinguished through  the  choice 


(5.1 5) 


Note  that  the  lines  of  principal  stress  are  not  uniquely  deter- 
mined in  the  case  of  stress  isotropy,  e.g.  when  regime  A or  C 
applies  throughout  a region. 

The  fundamental  equations  will  now  be  developed  in 
terms  of  the  geometry  of  the  a, p - co-ordinate  system  (see  Fig. 
6)r  Let  ih ( n.  A)  be  the  angle  made  by  the  tangent  to  the  a-line 
through  the  point  P(a,  {3)  with  a fixed  direction  5 and  let 
ha(a,f3)fta  and  h^(a,  f3)  6f3  be  the  elements  of  length  on  the 


4 r;.  , N \ ' 
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a-  and  p-lines,  respectively.  Let  x„  and  xQ  be  the  (algebraic) 

a p 

curvatures  of  the  a-  and  (3-lines,  and  be  defined  by 


V - -L.  - JL  Q _ 1 1 d± 

"«  “ Pa  ^ ta  ’ *P  " ^ ip-  (5-l6) 

Simple  geometry  applied  to  the  quasi-rectangular  element  PQRS  of 

Fig.  7 shows  that 

ah  3hn 


r~  = h h0x  , — & - h hftxQ . 
ap  a p a 9 0a  a P P 


(5.17) 


2 „ 2 

Further  since  -•  it  follows  that  the  curvatures  are 

related  through  the  condition 


axft  p 

+ c-§-  + x„ 


h^fe  + Ha  + *p 


Xr2  = 0. 


(5.18) 


Let  va  and  be  the  components  of  velocity  in  the 

directions  of  increasing  a and  p,  respectively.  Then  it  is 

easy  to  show  that  the  corresponding  strain-rates  e , gn,  v n 

<*  p7  'ap 

given  by 


- 9Va 
ea  h38a 


+ x„v 


ayp  ’ ep  = h^p  + Hpva 


3Vp  9va 

Yap=  V«  + hj^  ~ Vc  - «pvp  • 


(5.19) 


The  forces  acting  on  a small  element  h 6a  x hR6p  x t 

ct  p 

are  shown  in  Fig.  7.  The  differential  equations  of  equilibrium 


tiro  uo v wa 


taVa5  - VpVp  - 0 > %(h«V  - VfV.  = °-  <5-20) 

and  after  use  of  Eqs.(5.17)  these  equations  talre  the  form 


■J  . + M ‘ , 

... 


•i 
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6a 


a * uco 

- + = 0 > -gp  + hpxa^|3  ~aa)  = °*  (5.21) 

The  rate  of  plastic  work  per  unit  area  of  middle  surface  is 

W = t«Jaea  +<5^  ).  (5.22) 

The  yield  condition  is  given  by  replacing  and  dg 
by  d a and  dp  respectively,  i.e. 


“ " dA  > “ d. 


(5.23) 


'p  » x2  _ua  Jo  • 

In  an  isotropic  material  the  principal  axes  of  stress  and  plastic 
strain-rate  must  coincide,  i.e. 


t,(S=  0 


(5.2W 


throughout  the  plastic  region.  This  condition  imposes  a very 
severe  condition  on  the  velocity  field.  Apart  from  relation 
(5.24)  the  flow  rule  is  given  by  the  following  relations. 


Regime 

ea 

6P 

Conditions 

A 

-(l-q)X 

-qX 

da  = dp  ---0  ; £ ^ q ^ 1 

A B 

0 

—X 

°<  da  <do  » dp  = 0 

3 

qX 

-(l-q)X 

°a  =do»dp=°5°<q<l 

BC 

X 

0 

da  =dc  , 0 < dp  <d0 

C 

(l-q)X 

q 

da  = dp  = 0O  > o 1 q £ i- 

(5.25) 

Here  X(a,  B)  is  non-negative  and  q(a,  (3)  is  subject  to  the  stated 
inequalities,  but  otherwise  these  quantities  are  unrestricted  as 
yet. 


r . < _■  • - j 


' j i * s vf1'  r y*  - ** 

\ i*-  i?'1*  ■ " 
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The  equation  of  continuity  is 


1 i 

ha  9a  + hp  ap 


+ xpva+  *av{3 


t* 

t 


(5.26) 


The  set  of  six  equations  (5.18),  (5.21)  and  (5.23)- 
(5.25)  in  the  six  unknowns  da,  dp,  xa,  xp , vQ  and  Vp  form  the 
basis  for  the  determination  of  the  stress  and  velocity  fields 
in  the  plastic  region.  Note  that  the  first  four  of  these  equa- 
tions, viz.  Eqs.(5.l8)»  (5.21)  and  (5.23)  involve  only  the 
quantities  da,  dp,  xa  and  xp.  The  other  two  equations  involve 
the  velocities  and  are  homogeneous  in  these  quantities. 

In  the  rigid  (elastic)  region  the  equations  of  equi- 
librium remain  the  same,  but  now  the  strain-rates  all  vanish 
and  the  yield  condition  is  replaced  by  the  requirement  of 
compatibility 


-55)  + 4(Ep  -§f)  }<0n 


+ dp  ) = 0. 


(5.27) 


The  set  of  four  equations  (5.18),  (5.21)  and  (5.27)  form  the 
basis  for  the  determination  of  the  four  unknowns  a , dQ5  x„and 

(X  o vA 

Xp  throughout  the  rigid (elastic)  region. 

An  essential  preliminary  to  devising  methods  of  in- 
tegration of  the  plastic  equations  is  to  determine  their 
mathematical  nature.  The  rigid  ( elastic) equations  are  of  course 
elliptic  in  type.  Suppose  that  the  stress  and  velocity  fields 
are  known  throughout  the  plastic  region  bounded  by  a simple 
closed  curve  T (see  Fig.  8).  For  simplicity  it  is  assumed  that 
T has  everywhere  a continuously-turning  tangent,,  but  this 


condition  may  easily  be  relaxed  to  the  simpler  assumption  that 
P is  composed  of  a finite  number  of  arcs  with  the  above  property. 
Let  P be  a typical  point  of  I'  ; and  let  n and  s denote  local 
rectangular  Cartesian  co-ordinates  with  origin  P,  the  directions 
of  increasing  n and  s being  along  the  outwards -drawn  normal  and 
positive  tangent  to  T.  In  what  follows  we  shall  consider  the 
set  of  equations  found  by  replacing  x and  y by  n and  s , 
respectively,  throughout  Eqs.(5.2),  (5.9)  and  (5.10).  The  form 
of  the  physical  equations  is  of  course  unaltered  by  the  choice 
of  this  new  system  of  co-ordinates. 

Now  the  values  of  the  components  of  the  stress  and 
velocity  fields  (satisfying  the  plastic  equations)  are  known 
throughout  the  interior  of  F.  Hence  the  values  of  these 
quantities  together  with  their  partial  differential  coefficients 
of  all  orders  - the  data  being  supnosed  analytic  - are  all 
known  inside  T.  We  now  ask  if  the  known  data  just  inside  T, 
together  with  the  supposition  that  the  region  just  outside  V is 
plastic,  is  sufficient  to  determine  uniquely  the  corresponding 
data  just  outside  T,  If  this  is  so  then  the  the  known  solution 
inside  T may  be  continued  a small  way  into  the  region  outside 
T,  and  so  on.  If  this  is  not  so,  save  when  additional  data 
are  assigned,  then  T will  be  a characteristic.  The  characteristic 
condition  corresponds  of  course  to  the  situation  when  the  field 
quantities  together  with  their  normal  and  tangential  derivatives 
are  not  necessarily  all  continuous  across  T.  In  other  words 
different  analytic  solutions  of  the  plastic  equations  can  only 
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touch  along  a characteristic,  and  the  contact  may  be  of  any 
order , 

First  note  that  the  results  (5.13)  apply  to  T,  i.e. 

a , t ; and  v are  continuous  across  T , but  a and  v may  be 

n 7 “ns  n 7 s s 

discontinuous  across  P , The  values,  say  of  ds  * just  inside 

and  just  outside  T will  be  denoted  by  a " and  a + , respectively 

s S' 

and  the  discontinuity  measured  positively  outwards  across  T , 
i.e.  d + - <3  is  denoted  by  [<j  ]:  and  similarly  for  other 

quantities* 

a)  The  stress  equations. 

The  stress  equations  will  be  considered  first.  The 
yield  condition  (see  Eq.(5.9)) 

Cj<30  - dn)(jd0  - dg)  = (5.28) 

applies  on  both  sides  of  T , and  hence  as  xns  is  continuous 

[(jdD  - dn)(Jd0  - ds)]  = 0.  (5.29) 

The  analysis  now  proceeds  differently  according  as  the  same 
value  of  j does  or  does  not  apply  on  both  sides  of  T. 

Case  (i)  > ,1~  = j+(=  ,1).  Equation  (5.29)  shows  that 

(dn  - jdc)[dg]  = 0.  (5.30) 

Now  if  the  regime  be  either  A or  C,  so  that  there  is  a homo- 
geneous isotropic  stress  field,  then  nothing  more  remains  to  be 
said.  Attention  may  therefore  be  confined  to  regimes  AB,  B and 
BC.  Then  d is  continuous  if  dn  £ jd0  but  is  not  proven  contin- 
uous if 
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when  t =0  and  hence  d and  d are  principal  stresses, 
ns  ns 


Suppose  for  the  moment  that  there  are  strong  dis- 


continuities in  the  stress  field  so  that  Eq.(5.31)  applies. 


Therefore  if  the  regime  is  AB  or  BC  then  T must  be  an  a-  or  a 


P-line,  respectively;  and  if  the  regime  is  B then  T must  be 
either  an  a-  or  a p-line.  The  stresses  dn  and  t , together 


with  their  tangential  derivatives  and  ~~~  , are  continuous 

0*3  cs 


across  r.  The  equations  of  equilibrium  (see  Eqs.(5.2)) 


^n  . dxns  A Q'^ns  , 06  s 
an  + as  ~ 0 * an  + as  = 0 


(5.32) 


apply  on  both  sides  of  T,  and  prove  that  is  continuous. 
Further  from  Eas. (?«28)  and  (5.31) 


(as  - n0>  i2  ■ 


(5.33) 


In  regimes  AB,  B and  BC,  the  principal  stresses  are  unequal,  and 


hence  d ^jd  . Therefore 

b U 


= o. 


(5,34) 


The  quantities  -g—  , ^~and  are  not  proven  continuous, and  frem 


Eqs.(?.32)  their  discontinuities  are  restricted  only  by 


L an  + as  J “ 


(5.35) 


Results  governing  higher  order  discontinuities  are  easily  found. 


r ^ i. 
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Now  suppose  that  weak  discontinuities  at  most  occur 
in  the  stress  field.  The  stress  d must  now  be  prescribed 

o 

continuous  when  Eq,(5.31)  applies.  All  three  stresses  dn,ds 

8dn  0dc 

and  together  with  their  tangential  derivatives 

3tns 

and  — — , are  continuous  across  T . Equations  (5.32)  prove  that  the 

normal  derivatives  and  ^Tns  are  continuous  across  T . The 

on  on 

only  remaining  first  order  derivative  of  the  stresses  not  yet 
known  is  -^£+,  end  is  given  from  Eq.(5.28)  by 

<«n  - - 0.  (5.36) 


The  conclusions  previously  established  in  respect  of  T for  d - 

S 

Art 

discontinuities  therefore  also  apply  for  - discontinuities. 

Results  governing  higher  order  discontinuities  are  easily  found. 

If  now  in  addition  is  prescribed  continuous  when  Eq.(5.31) 

on 

applies,  then  the  lowest  order  partial  differential  coefficient  tfiich 


admits  discontinuities  is 


*7  , and  this  only  when  Eq.(5.31) 

3n" 


applies. 

Thus  it  appears  quite  generally  that  discontinuities 
can  occur  only  if  the  characteristic  condition  (5« 31)  is  satisfied. 
The  stress  equations  are  parabolic  for  regimes  AB  and  BC  with 
the  a-  and  p-lines,  respectively,  as  characteristics;  and  are 
hyperbolic  for  regime  B with  the  a-  and  (3-lines  as  characteristics. 

Case  (ii)  i ,1~  £ ,1*.  Equation  (5.29)  shows  that 


do(do-°n'<3s")  + (dn  “ 3+do)[ds]  = 0. 


(5.37) 


The  detailed  analysis  is  omitted,  and  only  the  essential  results 
will  be  stated.  Note  first  that  regimes  A and  G cannot  meet 
because  this  would  involve  the  discontinuity  of  both  principal 
stresses.  If  r is  not  in  regime  B then  a strong  discontinuity 
in  ds  of  known  amount  always  occurs  and  is  accompanied  by  weak 
discontinuities  in  the  stresses  also  of  known  amounts.  If  V is 
in  regime  B then  a strong  discontinuity  in  o and  weak  discontin- 
uities  in  the  stresses  may  occur. 

b)  The  velocity  equations. 

The  velocity  equations  will  now  be  considered.  Let  the 
two  orthogonal  families  of  lines  of  maximum  shear  stress  be 
identified  through  the  orthogonal  curvilinear  co-ordinates  y and 
6 . A member  of  one  family  is  then  defined  by  an  equation 
6 = const,  and  is  called  a y-line;  and  similarly  a member  of  the 
other  family  is  defined  by  an  equation  y = const,  and  is  called 
a 6-line.  Then  at  all  points  of  the  stress  field,  let  the 
directions  of  increasing  y and  6 be  the  internal  and  external 
bisectors  of  the  directions  of  increasing  a and  B, 

It  has  been  proved  that,  although  vn  is  continuous 
across  F,  vs  is  not  necessarily  continuous  across  T.  Now  such  a dis- 
continuity in  vg  when  properly  regarded  must  correspond  to  high 
shear  strain-rates  near  T . Hence  as  the  principal  axes  of 
stress  and  strain-rate  coincide,  r is  either  a y-  or  a 6-line. 

Suppose  for  the  moment  that  there  are  strong  discon- 
tinuities in  the  velocity  field  so  that  Tis  either  a y-  ora  6-line. 
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The  velocity  component  vn,  together  with  its  tangential  derivative 
dvn 

, are  continuous  across  r.  The  flew  rule  (see  Eqs.(5.9) 


and  (5.10)),  except  at  a singular  point,  is 


• &«>*,-„>  • }<&  * £■>*, 


(5.38) 


which  implies  that 


(5Z»  ♦ >)/<2!d 

an  3s  3n 


2t  /(d  - d ), 

ns  n s 


ft*  *\  r>  \ 

v j • j y i 


Equation  (5.39)  is  of  course  the  condition  that  the  principal 


axes  of  strain-rate  and  stress  coincide.  Now  in  the  present 


case 


d„  = d , 

n s’ 


(5.4o) 


and  hence  from  Eq.(5.39) 


avn  avs 


an  “ as 


(5.4i) 


The  continuity  of  d correctly  implies  that  of  o , The  law  of 

* 1 b 


mass  conservation  (sec  Eq.(5.11))  gives 


av  avs 

an  as  z ’ 


(5.4  2) 


and  hence  from  Eq.(5.4l) 


£ln  _ £ls  _ _ le 
an  as  2 z* 


(5. >+3) 


Therefore 


- Hez  1 » 


(5.44) 


- i'ir* , > / >Sk7:-£&V  -.-Cl.  * i-*. 
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av 


and  so 


n 


av. 


an  » as 


ana  e„,  are  all  cither  continuous  or  discontin- 
2 

9v„  avc  av^ 


uous.  Accordingly  the  quantities  — - £ , — £ and 


are  not 


proven  continuous,  and  their  discontinuities  are  restricted  only 

by  Eq. (5.44).  The  discussion  of  the  stress  equations  has  shown 

that  if  the  sane  value  of  j applies  on  both  sides  of  r.  then 

there  can  be  no  associated  discontinuities  in  the  stress  field. 

Otherwise  Eq.(5.37),  remembering  that  d is  continuous,  shows 

s 

that  dQ  = dn  + ds  = d + and  hence  T is  in  regime  B.  Weak 

ct  p 

discontinuities  may  now  occur  in  the  stress  field. 

Now  suppose  that  weak  discontinuities  at  most  occur  in 
the  velocity  field.  The  velocity  vs  must  now  be  prescribed 
continuous  when  r is  a Y“  or  &-line.  The  two  velocity  components 

0Vn 

vn  and  vs  , together  with  their  tangential  derivatives  and 

0Vg 

, are  now  continuous  across  V , and  it  follows  from  Eq.(5.42) 


0V. 


that  and  e_  are  either  both  continuous  or  discontinuous  across 
on  z 

T . More  nrecisely, 


0 v_ 

" - C'.J  • 


(5.45) 


The  only  remaining  first  order  partial  differential  coefficients 

av- 

of  the  velocity  components  not  yet  considered  is  - , and  in- 

0n  1 

formation  concerning  this  quantity  must  be  found  from  the  flow 
rule.  Now  the  flow  rule  (5.38)  may  be  re-written  in  the  form 


avn 

on 


= \ ( .1  d -d  ) 

~o  s 


av 

a T 


s _ 


- \(jd0-dn)  , |(g^2-  + g~)  - (5.46) 


ns 


l '' 
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where  X (non-zero)  is  an  undetermined  factor  of  proportionality. 

9v 

Since  Trns  and  are  continuous  across  r , the  third  of  Eqs.(5.46) 


9v, 


shows  that  and  X are  both  either  continuous  or  discontinuous 
across  P * More  precisely, 

avc 


C— S-1  = 2t  [X]  . 

uii  ns 


(5A7) 


Let  us  assume  for  the  moment  that  e is  continuous  across  P . 
av,  ov- 

Then  as  both  •— u.  and  *l  are  continuous,  the  first  and  second  of 
9 s an  ’ 

Eqs.(5.46)  show  that  ) ar.d  X(jo  -0  ) are  continuous. 


n 


Now  d is  continuous  and  so 
n 


d0l>J]  = C^SJ  = da[X] 


(5.48) 


The  quantities  j and  d are  not  proven  continuous.  However  if 

s 

j is  continuous  then  the  characteristic  condition  (5. 3D  applies, 

and  r is  an  a -era  P-line  depending  on  the  particular  plastic 

regime  involved;  and  if  d is  continuous  then  d^  = d and  P is 
7 s ns 

either  a v-  or  a 6-line.  On  the  other  hand  if  a is  discontinuous 

_ z 

9vn  3vq 

across  P , then  is  discontinuous  and  is  continuous  as 

dn  os 

before.  In  this  case 

[X(jd0-ds)]  - - [sz]  , [X(jd0-dn)]  = 0 (5.49) 

or,  alternatively, 

dQ[Xj]  = [Xds  - ezJ  = dn[X]  . (5.50) 

Again  the  quantities  j and  d are  not  proven  continuous.  However 
if  j is  continuous  then  the  characteristic  condition  (5. 3D  still 


Aii-IIn- 
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applies ; and  if  d is  continuous  then 

s 


<VVW  - tez]  - 


(5.51) 


and  no  further  statement  can  be  made.  It  has  therefore  been 
proved  that  strong  discontinuities  in  the  stress  field  may  be 
associated  with  weak  discontinuities  in  the  velocity  field. 

The  essential  results  are  therefore  as  follows. 

Strong  discontinuities  in  the  velocity  field  can  occur  only 
across  a v-  or  6 -line,  and  are  not  associated  with  strong  dis- 
continuities in  the  stress  field.  Weak  discontinuities  in  the 
velocity  field  may  or  may  not  be  associated  with  discontinuities 
in  the  stress  field.  In  particular  if  attention  is  confined  to 
just  one  plastic  regime  then  in  the  former  case  the  characteristics 
of  the  stress  and  velocity  equations  coincide  whereas  in  the 
latter  case  the  velocity  equations  are  hyperbolic  with  the  y- 
and  6-lines  as  characteristics. 

Now  that  the  existence  of  discontinuities  in  the 
velocity  field  has  been  established  it  is  more  straightforward 
to  approach  the  question  in  terms  of  the  geometry  of  the  a,p  - 
co-ordinate  system.  In  all  cases,  for  simplicity,  attention  is 
confined  to  the  case  when  the  same  plastic  regime  applies  on 
both  sides  of  P. 

The  following  results  now  follow  from  the  flow  rule 
(5.25),  and  govern  weak  discontinuities  in  the  velocity  field. 
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Regime 


r i I 

C*q-»  I 

-[(l-qJX]  > 


C^J 

m • 


-[qx] 

-fx] 

-CCl-q)X] 


C(l-q)X]  f C<1X] 


(5.52) 


1 


Here  it  is  supposed  that  [eal  and  [e^]  are  not  both  zero,  and 
the  results  for  the  three  cases  follow  immediately. 


[•«’ 


Regime 
A,  B,G 


Conditions 
'[X]/  [qX]  ^0 

[\yo,  [qX!  ='n~ 

[x]*  [qx]  i 0 
[x]#> 

[XJ=  [qx] 

[x]*> 

[X]*)  , [qX]  = 0 


(5.53) 


Thus  if  discontinuities  in  both  e_  and  eQ  occur  then  regimes  AB 

a p 

and  BC  are  excluded;  if  a discontinuity  occurs  only  in  ea  then 
regime  AB  is  excluded;  and  if  a discontinuity  occurs  only  in 
then  regime  BC  is  excluded.  The  conditions  under  which  such 
weak  discontinuities  in  the  velocity  field  may  occur  have  been 
discussed  previously. 


• --  -V.  - 

• L- v ■ .1?. • >•.  * " i-V'-L 


Now  consider  strong  discontinuities  in  the  velocity- 

field.  Here  it  is  supposed  that  [ v 1 and  [vD]  are  not  both 

a p 

zero,  and  hence  at  least  one  of  ea  and  becomes  very  large 
near  T.  Thus  X must  become  very  largo  near  F . Let 

T = lim  I \(n)dn  <>  0).  J = lim  f q(n)  X(n)dn  (0  < J < I)  (5. 5*0 

.1 


where  the  limit  is  taken  as  the  thickness  of  a small  region, 
enclosing  F and  across  which  the  velocity  changes  sharply,  is 
made  to  tend  to  zero.  Now  r is  either  a y-  or  a 6 -line,  and 
accordingly  the  results  are  as  follows  subject  to  a certain 
convention  to  be  described  below. 


Regime 

± Oa]//? 

± [va]//2 

r 

A 

-U-J) 

-j 

AB 

0 

-i 

B 

J 

-(i-j) 

BG 

I 

0 

C 

I-J 

j 

(5.55) 


If  the  direction  of  increasing  n lies  in  the  first,  second,  third 
or  fourth  gradrants  formed  by  the  a-  and  P-lines  then  the 
quantities  [ vff  ] //?  , [ v^ ] /V?  are  to  be  prefixed  by  both  positive, 
negative  and  positive,  both  negative,  and  positive  and  negative 
signs,  respectively.  The  results  for  the  three  cases  follow 
immediately. 
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[vl 

u J 

[Vq] 

Regime 

Conditions 

to 

to 

A,B,C 

I t J t 0 

to 

0 

A,  C 

I t 0,  J = 

B 

I = J / 0 

I 

BC 

1^0 

0 

^0 

A,  C 

I = J t o 

AB 

It  o 

I 

B 

I t o,  J - 

(5.56) 


Thus  if  discontinuities  in  both  vQ  and  Vp  occur  then  regimes  AB 
and  BC  are  excluded;  if  a discontinuity  occurs  only  in  va  then 
regime  AB  is  excluded;  and  if  a discontinuity  occurs  only  in  Vp 
then  regime  BG  is  excluded.  In  all  cases  the  integration  of  the 
equation  of  continuity  across  the  line  of  discontinuity  shows 
that 


i i 

•'  i 
I i 

r i 


<±  [va]  ± [vp3)  - -lim  j szdn, 


(5.57) 


the  sign  convention  being  as  before  and  the  limit  being  taken  in 

the  usual  way.  Thus  as  previously  noted,  e is  not  in  general 

z 

finite.  The  remarks  made  following  Eq.(4.l5)  are  again  applicable 
This  completes  the  discussion  cf  discontinuities  in 
the  velocity  equations. 

It  is  now  possible  to  take  up  the  point  in  the  analysis 

of  the  flange  equations  where  discontinuities  in  ^ and  ^ were 

seen  to  lead  to  discontinuities  in  (t  i,  and  (d„).  i.  , 

xy  y-$o  y.  y=?d  ’ 


- i **  • 

j . ' 


respectively  (see  p.  22  ).  Clearly  the  former  discontinuity  now 
appears  to  be  inadmissible  whereas  the  latter  discontinuity  is 
admissible  if  the  web  edge  is  an  a- cr  a ,6-line  depending  on  the 
plastic  regime.  Finally  it  is  necessary  to  consider  discontin- 
uities in  the  flange  axial  and  angular  velocities.  If  Up  is 
discontinuous  then  the  velocity  component  vx  in  the  web  must 
also  be  discontinuous.  This  can  only  be  accommodated  through 

an  infinity  in  the  transverse  strain-rate  e_  satisfying 

z 

[vx]  = - lim  j*  ezdx  , (5.58) 

dv™ 

the  limit  being  taken  in  the  usual  sense.  If  is  discontin- 

8v 

uous  then  the  quantity  — JL  in  the  web  must  also  be  discontinuous. 

8x 

Such  weak  velocity  discontinuities  in  the  web  have  already  been 
discussed.  In  particular  if  the  plastic  regime  is  continuous 
then  the  web  edge  is  an  a-craP-line  depending  on  the  plastic 
regime. 

The  general  problem  of  the  integration  of  the  web 
equations  subject  to  various  types  of  boundary  conditions  will 
not  be  examined  here.  However  this  Section  will  be  concluded 
by  showing  that  the  stress  and  velocity  fields  possess  certain 
very  striking  properties.  The*  may  be  compared  with  certain 
results  for  plastic  plane  strain  fields  (see  C 11  ]). 

Consider  first  the  regular  plastic  regimes. 

1)  Regime  AB:  0 < da  < > gB  ~ Q* 

The  first  of  Eqs.(5.20)  is  immediately  integrable,  and 
shows  that  h^d^  is  at  most  a function  cf  {3  only;  i.e.  hgda  is 


An  nL 

nxi  — • 


^8 


constant  along  an  a-line,  but  in  general  the  value  of  this 

constant  will  vary  from  one  a -line  to  another.  The  other 

equation  shows  that  xa  is  zero,  i.e.  the  a-lines  arc  straight. 

Therefore  kq  = l\h'(Q),  and  hence  p d is  constant  along  an 

P np  rp  a 

a-line.  The  geometry  of  the  stress  field  is  therefore 
particularly  simple,  and  is  shown  in  Fig.  9(a).  Conversely 
any  geometric  field  of  this  type  may  be  associated  with  a stress 
field.  The  flow  rule  equations  (5.24)  and  (5.25)  show  that  the 
a-lines  are  lines  of  zero  rate-of -extension  and  that 
av0 


V v?  + V«  * °>  tjjh  + idfe  - Vp 

Tho  a-lincs  are  said  to  form  a 'fan' . If  the  evolute  of  the 
a-  lines  degenerates  to  a point  then  the  fan  is  described  as 
'centered*.  In  this  case  choose  a and  p to  be  polar  co-ordinates 
r and  6 as  shown  in  Fig.  9(b).  Thus 


dv. 


3v 


= o.  (5.59) 


a = r,  p = 9,  hQ  = 1,  hp  = r,  *a  = 0,  = r ; 

and  the  stresses  and  velocities  are  given  by 


(5.6o) 


rdp  = f ( 9)  , dQ  = 0 , 

\ = g(e)  , vQ  = g* (9) 


—.V.  f Q'\ 
* **  \ v / 


(5.61) 


where  f,  g and  h are  arbitrary  functions  subject  only  to 


g"(9)  + g (9)  + rh' (9)  < 0. 


(5.62) 


To  sum  ups  the  a-lines  are  characteristics  for  the 
stress  equations,  and  are  straight  linos  with  zero  rate-of- 
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extension;  and  along  an  a-line  the  stress  dQ  is  proportional 
to  the  curvature  of  the  p -lines  where  they  cross  this  a-line. 

2)  Regime  BC:  dQ  = dQ,  0 < d&  < dQ. 

It  is  straightforward  to  obtain  results  similar  to 
those  above.  It  is  found  that  the  p-lines  are  straight,  and 
that  Pa (d0-<5p ) is  constant  along  a (3-line,  the  value  of  this 
constant  in  general  varying  from  one  line  to  another.  The 
geometry  of  the  stress  field  is  shown  in  Fig.  10(a).  The  flow 
rule  shows  that  the  p-lines  are  lines  of  zero  rate-of-extension, 
and  that 


av 


h±  + Vp  * 0 ’ t 


av 


a_  aVp  — v v_  = o (5.63) 

sr  vT^7  Ka  a w* 


‘ap  h^ap  ha3a 
The  p-lines  form  a fan,  and  if  this  degenerates  to  a centered  fan 
(sec  Fig.  10(b))  then  the  stresses  and  velocities  are  given  in 

terms  of  polar  co-ordinates  r and  6 by 


r(do-0r)  = f<©)  , c5q  = do 

vr  = g(0),  vQ  = g*(0)  + rh(0) , 


(5.64) 


where  f,  g and  h are  arbitrary  functions  subject  only  to 


g"(6)  + g(0)  + rh* (©)  > 0 


(5.65) 


To  sum  up:  the  p -lines  are  characteristics  for  the 

stress  equations,  and  are  straight  lines  with  zero  race-of- 
extension;  and  along  a p-line  the  modified  stress  dQ-dp  is 
proportional  to  the  curvature  of  the  a-lines  where  they  cross 
this  S-line. 


..--4  '% 
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Now  consider  the  singular  plastic  regimes. 

3)  Regime  As  = 6a=  0. 

2 e 

All  points  in  the  stress  field  are  isotropic,  and 
there  are  no  unique  lines  of  principal  stress. 

4)  Regime  Bs  6 = <5n  , oQ  = 0. 

a v. . .J2 

The  equilibrium  equations  (5.20)  show  that  and  x^ 
are  both  zero,  i.e.  the  a-  and  p-lines  are  straight.  The  lines 
of  principal  stress  therefore  form  a rectangular  mesh. 

5)  Regime  C:  = d^=  ft. 

All  points  in  the  stress  field  are  isotropic,  and  there 
are  no  unique  lines  of  principal  stress. 

6.  Wagner  Tension-Field  Beam  Theory. 

It  will  now  be  shown  that  Wagner's  tension-field  beam 
analysis  [2]  is  a special  case  of  the  analysis  developed  in  the 
previous  Sections.  For  simplicity  attention  will  be  confined  to 
the  beam  ABGD  of  Fig.  11  whose  function  is  to  transfer  shear  load 
from  the  rigid  section  A3  to  the  fixed  rigid  section  CD,  The 
flange  connections  at  A,R,C  and  D are  supposed  to  be  pin-jointed. 


Wagner ' s theory  determines  bno  limi  t,  load,  when  trie  design  is 
such  that  collapse  occurs  through  failure  only  of  the  web. 

Let  the  angular  velocity  of  both  flanges  be  oj  . 
Consider  the  following  velocity  field  for  the  web: 


vx  = 0 > vy  = ’ ux- 


(6.1) 


. b*  >1*1/  ‘ 
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Then 


e - e = 0 , y = - oj  , 
x y ’ T xy  * 


(6.2) 


and  hence  the  principal  axes  of  strain-rate  (and  stress)  make 
angles  of  45°  with  the  horizontal  and  vortical  directions.  Thus 


:a  = i u,  ep  = - i w,  Yap“ 


(6.3) 


and  as  u > 0 it  follows  that  the  web  is  in  regime  B,  i.c., 


da  “ dc  » = °* 


(6.4) 


In  other  words,  the  waves  run  at  45°  to  the  horizontal.  Further 
from  Eqs.(6.4)  it  follows  that 


°x  * °y  ’ - V = * ao  • 


(6.5) 


ine  equations  for  the  forces  and  bonding  moments  in  the  flanges 


— — + i td  — 0 dQ  + 1 +-H  — q dM  _ a 
dx  - ? too  u»  5x  ± 1 td°  “ U’  dx  " Q " °» 


(6.6) 


the  upper  and  lower  signs  corresponding  to  the  upper  and  lower 
flanges,  respectively.  The  horizontal  pull  by  the  web  on  the  end 
AB  is  dQtd,  and  this  must  be  reacted  by  compressive  forces 
id0td  in  each  flange.  Remembering  that  M(0)  = Mot)  = 0 for 
each  flange  it  then  follows  from  Eqs.(6.6)  that 


■ tip'1"  n t t) 
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F = + £ to0M-x)  - i dQtd,  I 
Q = + i tdoa-2x), 

M = + £ tdQ  xM-x). 


(6.7) 


The  downwards-acting  shear  force  applied  to  the  end  section  A3 


P = £ d0td. 


(6.8) 


The  velocity  field  (6.2)  requires  that  the  flanges  are 
rigid,  and  hence  the  yield  condition  for  the  flanges  is  not 
violated.  Now  M is  non-negative  and  non-positive  for  the  upper 
and  lower  flanges,  respectively.  Thus  it  is  necessary  that 


F2/Fd2  + K/M0  -lio, 


(6.9) 


i.e. 


t2{+(^— x)  - £d}2/l6a2b2  + tx(£-x)/4a2b  - 1 < 0.  (6.10) 


| 


The  maximum  value  of  the  left-hand-side  of  (6.10)  occurs  when 


x/t  = |l  - ^ U + J|)}/(2  - t/2b) . 


(6.11) 


Thus  if 


i1  * &C1 ' + ~ ?d4  /!1  >0  (6a2> 


then  the  beam  will  fail  in  shear  at  the  load  given  by  Eq.(6.8). 

An  important  extension  would  be  to  apply  the  present 
analysis  to  tapered  tension— field  beams. 


. . ' ■>../  '•1  i 


Under  a given  type  of  loading  a given  structure 
usually  fails  in  a well-defined  way.  A minimum  weight 
structure,  however,  is  capable  of  failing  in  any  one  of  a 
number  of  ways.  The  limits  of  economy  of  material  in  structures 
have  not  yet  been  investigated  in  a systematic  and  comprehensive 
manner.  Some  references  to  previous  work  have  been  given  by 
Hopkins  and  Prager  [23  ] . 

One  of  the  immediate  steps  required  in  the  development 
of  the  present  investigation  is  to  discuss  the  limits  of 
economy  of  material  in  a tension-field  beam.  Some  brief 
introductory  remarks  to  this  problem  will  now  be  given. 

The  present  discussion  will  be  limited  to  the  beam 
shown  in  Fig.  11.  This  beam  involves  three  structural  elements, 
viz.  the  upper  and  lower  flanges  and  the  web.  The  type  of 
failure  that  occurs  in  the  beam  depends  upon  the  relative 


Ac  , = compression  and  tension  flange  cross- 

sectional  areas. 
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Then 

G = flit  | (Ac  + At)/dt  + 1 

it  being  assumed  that  A and  A,  are  constant.  In  the  case 

c t 

discussed  in  Section  6 the  limit  shear  load  P is  ? 0Qdt,  and 
hence  Eq.(7.1)  can  be  written  in  the  form 

G = Y^1 * 3  | (Ac  + At)/dt  + * (7<2) 

More  generally  the  weight  may  be  written  in  the  form 

G = ay£3(_2|?))  (7.3) 


and  the  coefficient  a will  be  called  the  weight  factor.  The 
values  of  y,  P and  t being  fixed,  the  present  problem  is  to 
determine  the  values  of  A.  , Af  , d and  t which  render  the 
weight  factor  a minimum.  The  relations  (6.12)  impose  two  con- 
ditions. Other  conditions  must  be  found  through  the  study  of 
competi.vo  modes  of  collapse,  e.g.  those  in  which  either  the 
upper  or  the  lower  flange  remains  rigid  at  collapse  of  the 
beam.  The  general  procedure  may  be  expected  to  parallel  quite 
closely  the  analysis  given  in  Ref.  [ 23  ] , and  no  further  dis- 


cussion of  the  problem  is  attempted  here. 
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Fig.  5 Co-ordinates  on  discontinuity  surface 
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Fig.  6 Geometry  of  the  a,/5- co-ordinate  system 
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